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Abstract 


Cloud microphysical and optical properties are inferred from the bidirectional 
reflectances simulated for a single-layered cloud consisting of an external mixture of ice 
particles and liquid droplets. The reflectances are calculated with a rigorous discrete 
ordinates radiative transfer model and are functions of the cloud effective particle size, 
the cloud optical thickness, and the values of the ice fraction in the cloud (i.e., the ratio of 
ice water content to total water content). In the present light scattering and radiative 
transfer simulations, the ice fraction is assumed to be vertically homogeneous; the habit 
(shape) percentage as a function of ice particle size is consistent with that used for the 
Moderate Resolution Imaging Spectroradiometer (MODIS) operational (Collection 4 and 
earlier) cloud products; and the surface is assumed to be Lambertian with an albedo of 
0.03. Furthermore, error analyses pertaining to the inference of the effective particle sizes 
and optical thicknesses of mixed-phase clouds are performed. Errors are calculated with 
respect to the assumption of a cloud containing solely liquid or ice phase particles. The 
analyses suggest that the effective particle size inferred for a mixed-phase cloud can be 
underestimated (or overestimated) if pure liquid phase (or pure ice phase) is assumed for 
the cloud, whereas the corresponding cloud optical thickness can be overestimated (or 
underestimated). 
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1. Introduction 


Both ice and water phase particles may coexist in clouds when the temperature is 
between -40°C and 0°C [1], [2], Such mixed-phase clouds are important for several 
reasons. First, the difference in saturation vapor pressure between ice and water means 
that these clouds evolve differently from single-phase clouds. Additionally, mixed-phase 
clouds alter the radiative flux differently from single-phase clouds, thereby necessitating 
different treatment of clouds in both numerical weather forecasting and general 
circulation models (e.g., [3]). Furthermore, mixed-phase clouds can create hazardous 
conditions for aviation due to decreased visibility and aircraft icing [4], [5]. 

Mixed-phase clouds are frequently observed in the atmosphere [1], [6]. Recent 
research on mixed-phase clouds has focused on the analysis of ground-based or aircraft 
measurements [1], [6], the microphysical characteristics in cloud processes [7], [8], and 
modeling simulations [3], [9]. Flowever, current operational satellite-based retrievals of 
cloud microphysical and optical properties are based on the assumption of a cloud being 
composed entirely of either ice particles or water droplets [ 10]-[ 12]. This study 
investigates the cloud optical and microphysical properties for clouds containing a 
mixture of both ice and water particles and compares them with retrievals that incorrectly 
assume that the phase is solely ice or liquid water. 

2. Radiative Transfer Model 

The radiative transfer model used in this study discretizes the atmosphere into 100 
layers. The cloud-free atmospheric optical thickness for each layer is computed using the 
MODTRAN radiative transfer model [13]. The bidirectional reflectance at the top of the 
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atmosphere is computed using the discrete-ordinate radiative transfer model (DISORT) 
[14] with a precalculated clear-sky atmospheric gaseous optical thickness vertical profile 
from MODTRAN. For a cloud layer, the bulk single-scattering properties of ice particles 
and water droplets are available from a precalculated library. The single-scattering 
properties of ice particles, covering a range of sizes from 1 to 4000 pm, are computed 
with the approach reported in [15]. The mean (or bulk) scattering properties are 
calculated using 30 in situ size distributions measured for a variety of midlatitude and 
tropical cirrus cloud systems obtained during several field campaigns [16], [17]. 

In this study, the particle habit (or shape) percentage as a function of particle size for 
a given size distribution of ice crystals is consistent with that used in the MODIS 
(Moderate Resolution Imaging Spectroradiometer) operational (Collection 4 and earlier) 
cloud products [18]. Specifically, the habits of ice particles considered for a given size 
distribution includes bullet rosettes, hollow columns, aggregates, and hexagonal plates. 
For small ice particles (D max < 70pm, where D max is maximum dimension), the ice habit 
mixture is composed of 50% bullet rosettes, 25% hexagonal plates, and 25% hollow 
columns. For large ice particles (D max > 70pm), the mixture is composed of 30% 
aggregates, 20% bullet rosettes, 20% hexagonal plates, and 20% hollow columns, 
following [19] and [12]. The single-scattering properties of water clouds are calculated 
based on the Lorenz-Mie scattering, assuming 50 different Gamma size distributions. 

The single-scattering properties of mixed-phase clouds are calculated following [20]. 
Here, the mixed-phase clouds are assumed to be a uniform mixture of ice and water 
particles. The total water content ( TWC ) is the sum of ice water content ( IWC) and liquid 
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water content ( LWC ), i.e., TWC = IWC + L WC. Following [20], the ice mass fraction (y), 
is defined as follows: 


IWC 
TWC ‘ 


( 1 ) 


The effective particle size D em for a mixed-phase cloud is derived in [20] by: 


H em 
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where D e and p indicate the effective size and mass density, respectively; and the 
subscripts i, w, and m in Eq. (2) denote ice, water, and mixed-phase clouds, respectively. 
The optical thickness of a mixed-phase cloud (t,„) is given by 

Tn =T,.+T W = A zf QeiMdD + A zf Q ew A w n w dD (3) 


where A z is the geometric thickness of the cloud, Q e is the extinction efficiency, A is the 
projected area of the particle, n is the size distribution of particle number density, D is the 
characteristic length of the particle. It can be shown that the optical thicknesses pertaining 
to the ice and water components of a mixed-phase cloud are given, respectively, by the 
following expressions: 


( Qei ')^emPmY 

(Q em )D eiPi Xm ’ 
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The detailed mathematical expressions for the quantities involved in Eqs. (4) and (5) can 
be found in [20], 

3. Methodology 
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We calculate the bidirectional reflectance at the top of the atmosphere using the 
radiative transfer model discussed in the previous section. The bidirectional reflectance is 
defined as 


n 0 ,(p 0 ) 


jeT(O;^,0) 


( 8 ) 


where 7(O;/i,0) is the reflected radiance at the top of atmosphere, Fq is the incident solar 
flux at the top of the atmosphere, po is cosine of the solar zenith angle, (po is the solar 
azimuthal angle, p is cosine of the sensor zenith angle, and (p is the sensor azimuthal 
angle. In this study, we simulate the bidirectional reflectance at high latitudes where 
mixed-phase clouds are frequently observed and use 55.8°, 7.0°, and 20.0° for the solar 
zenith angle, the sensor zenith angle, and the difference between solar and sensor 
azimuthal angles (/. e., the relative azimuth angle), respectively. The corresponding 
scattering angle is approximately 130°. 

A popular method pioneered by [21], [10], and [22] infers cloud microphysical and 
optical properties by using the bidirectional reflectance observations in two or more 
spectral bands that have varying amounts of cloud particle absorption, e.g., a 
nonabsorbing visible band and an absorbing near infrared band, or even two shortwave 
infrared bands [11]. The reflectance in a nonabsorbing band (e.g., 0.67, 0.87pm) is 
largely dependent on cloud optical thickness, whereas the reflectance in an absorbing 
band (e.g., 1.6, 2.1, and 3.75 pm) is sensitive to the cloud particle size [23], In this study, 
we use a nonabsorbing and absorbing band defined by the MODIS spectral bandpasses at 
0.87 and 2.13 pm, respectively, which is a combination used for the primary MODIS 
operational cloud retrieval [18]. 
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This method is most accurate over a dark ocean surface away from sunglint, where 
the reflection comes from the cloud layer with little influence from the underlying surface 
[12]. On a bright surface such as ice or snow where the surface reflectance is large, the 
signal from the cloud is mixed with the signal from surface. In this study, a surface 
albedo of 0.03 is used both at 0.87pm and 2.13pm [24], [12] which is a value for dark 
ocean surface. Sunglint is not a factor for the solar-viewing geometry defined in this 
study. 

4. Results 

Figure 1 shows the variation of the effective particle diameter (Fig. 1(a)) and optical 
thickness (Fig. 1(b)) for clouds composed solely of ice, water, and a range of ice mass 
fractions. The water cloud has an effective particle diameter, D ew , of 15pm and an optical 
thickness of t w = 10. The ice cloud has an effective particle diameter, D ei , of 50 pm and an 
optical thickness of t; =10. In Fig. 1(a) when the ice fraction y=(), i. e . , a pure liquid water 
cloud case, the effective particle diameter becomes that of the water cloud. When y=l, i. e . , 
a pure ice cloud case, the effective diameter becomes that of the ice cloud. For values of y 
between 0 and 1, the effective particle size varies smoothly between the ice and water 
cloud value. When y=0.5, for example, the ice water content equals the liquid water 
content, but the effective particle diameter (23.5 pm, see Eq. 2) is less than the arithmetic 
mean of the ice and water cloud effective particle sizes (32.5 pm). This occurs because 
there are larger numbers of smaller particles in the water cloud. 


6 



Fig. 1(b) shows similar results as for Fig. 1(a), but for a cloud with total optical 
thickness, r m , equal to 10. Note that the magnitude of the derivative of the r w and t, 
components increases rapidly with y. 

Figure 2 shows the response of the bidirectional reflectance at 0.87 and 2.13pm to 
changes in y. These results are superimposed on theoretical isolines of r and D e for pure 
ice (Fig. 2(a)) and pure water (Fig. 2(b)) clouds. Each symbol (+) corresponds to the ice 
fraction, which changes from 0 to 1 in increments of 0.1. The reflectances change 
gradually between values corresponding to a homogeneous water and ice cloud. The 
reflectance at 2.13pm for a water cloud (y=0) is larger than that for the ice cloud (y=l) 
because (a) ice is more absorbing than liquid water, and (b) ice cloud particles tend to be 
larger than for the water cloud and therefore have a smaller cloud reflectance. Note that 
the cloud reflectance tends to decrease with increasing particle size. Thus, as y increases, 
the effective particle size inferred from the liquid water theoretical curves increases. The 
optical thickness in Fig. 2b also increases with y because the ice particles have a larger 
asymmetry parameter due to their larger size and nonspherical habits. 

Figure 3 compares the differences between homogeneous (with respect to 
thermodynamic phase) ice and water clouds for effective particle diameters (Fig. 3(a)) 
and optical thicknesses (Fig. 3(b)) retrieved from the theoretical relationships of Fig. 2. In 
this example (D ( . ; =50pm, D ew =\ 5pm), the effective particle diameter inferred with the 
assumption of a homogeneous water cloud is smaller than the actual particle size of the 
mixed-phase cloud. Conversely, the effective particle diameter obtained assuming an ice 
cloud is larger than that of the mixed-phase cloud. Figure 3(b) shows the differences in 
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the inferred optical thickness. For this case (r m =10), the optical thickness obtained with 
the water (ice) cloud assumption is larger (smaller) than that of the mixed-phase cloud. 

Figure 4 shows the variation of the bidirectional reflectance at 0.87 and 2.13pm with 
both y and optical thickness for mixed-phase clouds as well as for homogeneous ice and 
water clouds. The effective particle diameter of the ice and water clouds is fixed at 50pm 
and 15pm, respectively. The optical thickness varies from 1 to 15. In Fig. 4(a) when y=l, 
the reflectances are those values corresponding to the optical thickness of 1, 3, 5, 8, and 
15 for ice clouds, respectively. In Fig. 4(b) when y=0, the reflectances are those values 
corresponding to the optical thickness of 1, 3, 5, 8, and 15 for water clouds, respectively. 
It is clear from this plot that substantial errors in inferred effective particle size and 
optical thickness parameters may be incurred if a mixed-phase cloud is treated as a 
single-phase ice or water cloud when using a bispectral retrieval algorithm. 

5. Summary 

Mixed-phase clouds are a frequent phenomenon in the atmosphere and are important 
for several reasons. In this study, we have investigated the accuracy of the retrieval of 
two properties, effective particle size and optical thickness, for a single layered cloud 
containing both ice and water particles using the bidirectional reflectance at two 
wavelengths. The mixed-phase cloud is modeled as an external mixture of ice and water 
particles rather than as two distinct layers, each of which has hydrometeors of a single 
phase. Bidirectional reflectances are computed using a discrete ordinates radiative 
transfer model over an ocean surface having an albedo of 0.3. Based on these simulated 
reflectances for both single-phase and mixed-phase clouds, the effective particle diameter 
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and optical thickness are inferred. The sensitivity of these inferred values to the ice 
fraction is investigated. It is shown that the assumption of a single thermodynamic phase 
can induce errors in the retrieved parameters when the cloud is composed of both ice and 
water, since the bidirectional reflectances are influenced by the ice fraction. Specifically, 
the effective particle size retrieved for a mixed-phase cloud from the visible and near- 
infrared bands can be underestimated (or overestimated) if pure water (or ice) phase is 
assumed, whereas the corresponding optical thickness of the mixed-phase cloud can be 
overestimated (or underestimated). While we present results for only one solar-satellite- 
viewing geometry, the sensitivity to ice fraction is expected to hold for viewing 
conditions that are unaffected by sunglint. Finally, we note that the identification of 
mixed-phase clouds as well as their discrimination from well-separated multiphase cloud 
layers (such as thin cirrus overlying a lower-level water cloud) remains a very 
challenging task for passive satellite sensors. 
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Figure legends 

Figure 1. The variation of (a) effective particle diameter and (b) optical thickness as a 
function of ice fraction (y) in a mixed-phase cloud. Subscripts i, w, and m refer to an ice, 
liquid water, and mixed-phase cloud, respectively. 

Figure 2. The variation of the bidirectional reflectance of mixed-phase cloud at 0.87 and 
2.13pm with the ice fraction (y) superimposed on the theoretical r-D e curves for the (a) 
ice cloud and (b) water cloud, y varies from 0 to 1 with the (+) symbols from top to 
bottom. 

Figure 3. The comparison of the (a) effective diameter and (b) optical thickness retrieved 
assuming a homogeneous ice and water cloud as a function of ice mass fraction (y). The 
true mixed-phase cloud parameters (r m , D em ) are also shown. 

Figure 4. Same as Fig. 2 except that the optical thickness of mixed-phase cloud varies 
from 1 to 15. 
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Figure 1 . The variation of (a) effective particle diameter and (b) optical thickness as a 
function of ice fraction (y) in a mixed-phase cloud. Subscripts i, w, and m refer to an ice, 
liquid water, and mixed-phase cloud, respectively. 
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Figure 2. The variation of the bidirectional reflectance of mixed-phase cloud at 0.87 and 
2.13pm with the ice fraction (y) superimposed on the theoretical r-D e curves for the (a) 
ice cloud and (b) water cloud, y varies from 0 to 1 with the (+) symbols from top to 
bottom. 
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Figure 3. The comparison of the (a) effective diameter and (b) optical thickness retrieved 
assuming a homogeneous ice and water cloud as a function of ice mass fraction (y). The 
true mixed-phase cloud parameters (r m , D em ) are also shown. 
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Figure 4. Same as Fig. 2 except that the optical thickness of mixed-phase cloud varies 
from 1 to 15. 
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